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Stability analysis of 5D gravitational solutions with N bulk scalar fields 
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We study the stability of 5D gravitational solutions containing an arbitrary number of scalar 
fields. A closed set of equations is derived which governs the background and perturbations of 
A'^ scalar fields and the metric, for arbitrary bulk and boundary scalar potentials. In particular 
the effect of the energy-momentum tensor of the scalar fields on the geometry is fully taken into 
account, together with all the perturbations of the system. The equations are explicitly written as 
an eigenvalue problem, which can be readily solved to determine the stability of the system and 
obtain the properties of the fluctuations, such as masses and couplings. As an example, we study a 
dynamical soft-wall model with two bulk scalar fields used to model the hadron spectrum of QCD 
and the Higgs sector of electroweak physics. It is shown that there are no tachyonic modes, and 
that there is a (radion) mode whose mass is suppressed by a large logarithm compared to that of 
the other Kaluza-Klein modes. 



PACS numbers: 



I. INTRODUCTION 



Brane-world models have proved to be very useful for 
model building, providing a way to not only address 
the hierarchy problem in the Standard Model [iHSI, but 
also to explain the hierarchy of fermion masses and mix- 
ings 043 ■ ^ variety of constructions are usually em- 
ployed with fields, other than gravity, in the bulk and/or 
on the branes. Furthermore, aided by the AdS/CFT cor- 
respondence, a dual description of strongly coupled four- 
dimensional (4D) gauge theories can be obtained 
These fields are normally treated as test fields, where the 
corresponding energy density is sufficiently small not to 
disturb the background geometry. Therefore, for a given 
geometry one performs a Kaluza-Klein (KK) decomposi- 
tion. Solving the field equations in this geometry gives 
the bulk profile and the eigenmasses of the KK modes. 
The profile, in turn, determines the couplings of the KK 
modes. 

Eventually, a complete model has to provide the dy- 
namical elements that determine the bulk geometry, and 
the brane positions (when present). As an example, the 
Randall-Sundrum model Q, does not solve the hierar- 
chy problem until the exact location of the TeV brane is 
fully specified. Most of stabilization mechanisms, start- 
ing from [TT| . typically employ one scalar field, due to its 
simplicity. 

The system of perturbations in the presence of a bulk 
scalar was studied, for example, in Ref. (see also 
[isj). This analysis considered the properties of various 



excitations coupled to SM fields, which were assumed to 
be localized on the brane. Of particular relevance was 
the study of the radion, which is the lightest of these 
scalar perturbations. The analysis of fT3| was limited to 
a regime of small backreaction of the bulk scalars on the 
bulk geometry. Later, generalized this study to a 
bulk field with arbitrary bulk/brane potential. 

However, there are systems in which more than one 
dynamical scalar field is relevant. For instance, in the 
soft- wall model of two bulk scalar fields are used to 
obtain a solution of the Einstein equations. Furthermore 
in string theory, multiple scalar fields (such as the dila- 
ton and moduli fields) are quite common. This motivates 
generalizing the study of [ij] to the case of N bulk (real) 
scalar fields, which will be presented in this paper. The 
generalization to N bulk scalar fields has also been pre- 
viously considered in Refs. [l6l - [l8j . In our analysis we 
make no assumption on the form of the bulk and bound- 
ary potentials. 

The perturbations are obtained by studying the lin- 
earized problem. It is useful to distinguish them as 
scalar/vector/tensor modes under Lorentz transforma- 
tions of the ordinary 4D spacetime: we assume that the 
noncompact space has a Minkowski metric. As a con- 
sequence the scalar/ vector/tensor modes are decoupled 
from each other at the linearized level and therefore can 
be studied separately. The mode count proceeds as fol- 
lows. There are 15 perturbations in the five-dimensional 
(5D) metric, and N perturbations from the bulk (real) 
scalar fields. Of these, 5 are nondynamical and 5 more 
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can be removed by a gauge choice^. This leaves 5 + N 
modes of which, 5 are in the tensor sector (transverse 
and traceless /i^y) and A'^ remain as scalar modes. There 
are no vector modes. Finally, there is a scalar mode for 
each brane (representing the fluctuation of its bulk po- 
sition); however, we show that these modes are not ex- 
cited as long as the branes have only a tension (namely, 
no other brane source with a different equation of state 
is dynamically relevant; this is the usual assumption in 
brane compactifications) and so we will ignore them. 

It is technically challenging to identify the N scalar 
degrees of freedom because they can arise either from 
the components of the metric perturbations or the ac- 
tual bulk scalar field perturbations. The number can 
be reduced to A'' by using constraint equations to derive 
a manageable set of closed equations for the N modes. 
This must be done for both the bulk and brane equa- 
tions. The latter equations are boundary conditions and 
we will consider both cases either with or without branes 
at the boundaries (in the second case, the boundary con- 
dition is typically given by the requirement of normaliz- 
ability). We derive this system of equations, which for 
convenience is written as an eigenvalue problem, with the 
appropriate number of equations needed to obtain a so- 
lution. Our goal is to provide an explicit formulation of 
the eigenvalue problem that can be used to study any 
5D model with an arbitrary number of bulk scalar fields. 
By solving the eigenvalue problem the physical proper- 
ties of the perturbations (masses and couplings) can be 
obtained. This can then be used to check the stability of 
the model, so that if < for some scalar modes, the 
background solution is unstable. 

While the eigenvalue problem may be solved ana- 
lytically for the simplest cases, in general a numerical 
method is needed to obtain a solution. We will employ 
the shooting method since the boundary equations at one 
boundary leave freedom for the choice of the mode func- 
tions, and of the mass eigenvalue. We show that, for TV 
scalars, the boundary conditions at one of the two bound- 
aries leave N unspecified quantities. The bulk equations 
are then used to evolve the solution to the other bound- 
ary, where there are precisely N constraints that must be 
satisfied (given by the boundary condition at this bound- 
ary). If these constraints are satisfied, then a physical 
mode has been found. 

As an example of the numerical method, we study the 
perturbation properties, and the stability problem, for 
the dynamical soft wall model of Ref. . The model is 
characterized by two scalar fields with a coupled poten- 
tial term. It is interesting because it leads to a KK mass 
spectrum with linear Regge-like trajectories, similar to 
the hadron spectrum in QCD. We conduct a dense scan 



in parameter space and find approximately 100 modes 
with positive m^, and no mode with negative rri^. More- 
over, we find an interesting behavior of the radion mass. 
The mass is suppressed by a large logarithm, that in the 
dual CFT interpretation corresponds to how scale invari- 
ance is broken by quantum corrections. This is the same 
suppression present in the Goldberger-Wise mechanism 
fll, 12]. 

The paper is organized as follows. In Section |ll] we 
introduce the class of models that we are studying; we 
present the Lagrangian, the background solutions, and 
the most general set of perturbations. In the following 
three Sections we study the scalar, vector, and tensor per- 
turbations, respectively. We identify the physical modes, 
and perform their Kaluza-Klein decomposition. We pro- 
vide the explicit closed set of equations which can be 
solved to obtain the mass and the bulk profile of these 
modes through a boundary value problem. In Section 
IVII we formalize this boundary value problem for the two 
sectors (scalar and tensor) that contain physical pertur- 
bations. In Section |yill we study the perturbations of 
the model in Ref. [l5| as an example of how to use our 
formalism to solve the eigenvalue problem. In Section 
IVIIII we briefly summarize our main findings. Some more 
technical steps are then given in the three Appendices. 



II. 5D MODELS 

The goal of this work is to provide the tools for study- 
ing the perturbations, and the stability, of a wide class of 
models with one extra dimension and A^ bulk scalar fields 
{(pi}. Specifically, we consider models characterized by 
the action 
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]\r3 1 



^ Gauge invariant variables can also be used. We choose the ex- 
plicit gauge choice here, for brevity of exposition, but the two 
approaches are completely equivalent. 



- E / d'xV^{2M^K]j + Ui^,)} , (1) 

branes 

where M is the 5D fundamental scale. The first line in 
contains the bulk terms, which are restricted to lie 
between two boundary branes, if they are both present, 
or else either between a boundary brane and infinity, or 
between ±cxi. The overall factor of 2 is adopted from 
14| . where the bulk was assumed to be Z2 symmetric 
across each brane, and the symmetry was accounted for 
by restricting the bulk integral only to one side of each 
brane. This notation will also be used here, even in the 
cases where one or both branes are absent; it is trivial to 
reabsorb this factor away by a rescaling of M, ipi, and 
V. We assume that the scalars have a standard kinetic 
term in the bulk. Notice that the choice of sign for the 
kinetic term corresponds to 77^,^ = diag(— 1, 1, 1, 1). If 
present, a bulk cosmological constant is included as a 
constant term in V. The second line in ([1} is instead the 
brane action; 7 denotes the induced metric on the brane, 
and the notation [K]j with subscript J denotes the jump 
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of the quantity inside the square parenthesis across the 
brane, which in this case is the extrinsic curvature K. U 
denotes the potential of the scalars on the brane (which is 
a function of the value that the ipi fields have at the brane 
location), and, if present, a brane tension is included as 
a constant term in U (note that when both boundary 
terms are present, we do not require them to be equal). 
In addition we disregard the possibility of brane kinetic 
terms for the scalar fields. 

From the action ^ we obtain the Einstein equations 
in the bulk 



Gab ~ 



Tab 

M3" 



(2) 



where capital Latin indices run over all the coordinates 
and the energy momentum tensor is given by 



Tab = dAfidB(Pi + qab 



We also obtain the N bulk equations for the scalars 
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dV 



dipi 



(3) 



(4) 



Note that not all of the equations in ([2]) and (|4]) are in- 
dependent, due to the Bianchi identities. 

When a boundary brane is present, extremizing the 
action ([T]) leads to the boundary conditions for the scalar 
fields 



dU 

as well as the so-called Israel conditions 



'M3 



(5) 



(6) 



where K^^ = K^^ — K 7^^ , and S^i, denotes the stress 
energy tensor on the brane 



(7) 



Note that Greek indices run over the usual 3-1-1 dimen- 
sions only. The explicit definition and computation of 
the induced metric and the extrinsic curvature are given 
in Appendix A. 



A. Background 

At the background level, we assume a factorizable ge- 
ometry with 4D Minkowski slices: 



(8) 



It follows that the background bulk scalars can have a 
nontrivial dependence on the extra coordinate only: 



(fl = (f>i (z) 



(9) 



The only nontrivial equations in ([2]) then arise from the 
diagonal fiv, and z components. We write here one linear 
combination of these two equations, and the zz compo- 
nent: 
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A^V , 



(10) 



where prime (') denotes differentiation with respect to z. 
The scalar equations (jlj give instead 



A 



b',-A^V,^0 , 



(11) 



where Vj = . It is easy to check that the first equation 
of pUj) can be derived by combining (|lll) and the second 
equation of (I10[) (this redundancy is a consequence of a 
nontrivial component of the Bianchi identity). 

From the background expressions of the induced met- 
ric and the extrinsic curvature given in Appendix A, the 
Israel conditions (O have only a nontrivial part propor- 
tional to T]^^: 



A' 
A^ 



U 



(12) 



while the boundary conditions ^ can be rewritten as 

(13) 



A 2 



where U^i = The upper (lower) sign at the right 

hand side of these two equations refers to a brane at the 
left (right) of the bulk interval. If one or both branes 
are absent, equations (fT2|) and (fT3|) can be replaced by 
different boundary conditions at spatial infinity along the 
bulk. 

In the study of the perturbations, we often make use 
of the background equations written in this Subsection 
in order to simplify the linearized equations for the per- 
turbations, without writing this explicitly each time. 



B. Perturbations 

It is convenient to characterize the perturbations ac- 
cording to how they transform under 4D Lorentz trans- 
formations. Due to the background symmetry, modes 
that transform differently under these transformations 
are decoupled at the linearized level, and can be studied 
separately in our analysis. We therefore have the follow- 
ing decomposition 



5^. = (z) [7;^, (1 + 2*) + 2d^d,E 

= A^{z)[d^B + B^] , 
555 = A2(z)[l + 2$] . 



(14) 
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The modes ^, E, B, $ are scalar (with respect to 4D 
Lorentz transformations); the modes E^, are vector 
(we impose that they are transverse, = = 0), 

and hfj^y is a tensor (imposed to be symmetric, transverse, 
and traceless, d^^h^^y = /i^ = 0). There are also addi- 
tional scalar modes: the perturbations, 5(j)i of the bulk 
scalars Lpi = (f>i{z) +54'i, and the perturbations, Q of the 

brane positions Zj-thbrane = ^background position + C,]{x^^) 

{j runs over the number of branes). All the perturbations 
are functions of both x'^ and z, except for Q which are 
functions of x^^ only. Note that the decomposition pi)) 
becomes ambiguous for massless scalar KK modes, but 
can be studied using the light cone decomposition of [l^ . 
We assume that there are no massless scalar modes in the 
cases of interest. 

We need to fix the freedom of general coordinate trans- 
formations. Under the infinitesimal transformation 

(with df^£,'^ = 0) the metric changes as 

SgAB ^ SgAB - gTB,c(^ - 9^0^% - &A • (16) 

We can use this relation to see how the various modes 
in (jl4p transform. In particular we obtain: 



E^ E-i , B ^ B 
^ Ef^ - , 



(17) 



where we have only given the transformations relevant for 
the present discussion. This allows us to set = in 
(1141) and removes the freedom of the transformations in 
(ITSI) characterized by Similarly we can also set i? = 
in (HH), which fixes ^, and choosing _B = then fixes 
Therefore we see that one can always choose the gauge 
E^ — E — B ^ 0; this completely fixes the freedom of 
the coordinate transformations ([T5|) . This leaves the sets 
of scalar j^*, Scpi, Q}, vector {i?^} and tensor {hfj,^} 
modes; these three systems are decoupled from each other 
at the linearized level, and we will study them separately 
in the following three sections. 



the index j on the displacement Q ] 
Appendix A, we obtain 



Using the results in 



[{A' (1 - $ + 2*) + A^'} Tj^y - AC^,] 



U 

M3 



A"" (1 + 2*) r/^. 



(18) 



(Note that U also has a perturbation part). This equa- 
tion has two tensorial structures that need to vanish inde- 
pendently. In particular, we find that [C,/ii/]j — 0. Using 
the Z2 symmetry at the brane, and the fact that C, is odd, 
this equation in turns gives 



□C = o 



(19) 



where □ denotes the d'Alembertian operator in 4D 
Minkowski space. As is well known for the single scalar 
case, the brane bending mode is not sourced by the scalar 
fields, and is just a decoupled massless mode in the cur- 
rent context. Therefore, in the following we simply dis- 
regard this brane mode(s). 

We are left with the TV 4- 2 scalar perturbations 
{5", $, (5(/)i}. As we discussed in the Introduction, there 
are only N physical perturbations in the scalar sector. 
In Subsection IIII Al we present the linearized bulk equa- 
tions for the scalars, and show how the two nondynamical 
modes can be eliminated from two constraint equations. 
We actually define N scalar combinations Vi, that cor- 
respond to the canonical variables of the system. In the 
following three subsections we then compute the bound- 
ary conditions for these variables for different relevant 
cases. 



A. Bulk equations 

We start by considering the scalar projection of the 
linearized Einstein equations ([2]) in the bulk. The nv 
components read 



(...)ry^,-($ + 2vE')^^^ = 
The off-diagonal part requires that 



(20) 



III. SCALAR PERTURBATIONS 

In this Section we write the linearized equations for 
the scalar perturbations. The main goal is to remove 
the nondynamical degrees of freedom. As we wrote in 
the last paragraph of the previous Section, we start from 
the system of perturbations j^*, $, (50^, Cj}: where i runs 
over the number of bulk scalars (N)^ while j runs over 
the number of branes. 

We can immediately show that the brane displace- 
ments are decoupled, and do not introduce any instabil- 
ity. They only enter in the boundary conditions. Specif- 
ically, let us assume that a brane is present and con- 
sider the linearization of ([5]) at that location (omitting 



^ = . 

2 



(21) 



From now on, we will enforce this constraint to eliminate 
The /i5 linearized equations then give 



J $ 

2 A 3Af3 



(22) 



Eqs. ([?T|) and ([H]) arc the two constraints that allow the 
remaining two nonphysical modes in the scalar sector to 
be eliminated after the gauge fixing. 
Next, we introduce the N combinations 



(23) 
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which generahze the canonical variable introduced in 
for the case of a single scalar field. These modes are the 
canonical variables of the system. After the conditions 
PT|) and are used to eliminate the nondynamical 
modes, one can show (see Appendix B) that the remain- 
ing bulk equations are equivalent to the following system 
of N equations for the N modes Vi: 



2A3 




'kVk 



3M3A'2 



2M3 



(24) 



This is the explicit set of equations for the dynamical 
scalar modes of the system. 



The final result is 
3A' 



branc 



2 A " \iM^A' 



A 



3A/3 
V2 



X (2720^ 



V2A' 




(27) 



which are indeed equivalent to (|26p (we show this in Ap- 
pendix C). We stress that these equations are valid at 
the brane location; brane quantities are evaluated imme- 
diately next to the brane, and whenever an upper /lower 
sign appears on the right hand side, it refers to a brane at 
the left/right of the bulk interval, respectively. Finally, 
let us clarify the role of the □ operator on the right hand 
side. The bulk equations p4)) allow for a factorizable 
solution 



Vi (a 



(28) 



B. Boundary conditions for finite brane potential 



with 



After disregarding C, only the term proportional to 77^^ 
survives in the Israel conditions p^ . Using (|2ip and the 
Z2 symmetry, we rewrite it as 



24 — $ A 

A2 r ^ 2A' 



(25) 



where the left hand side is evaluated in the bulk imme- 
diately next to the brane, and the upper (lower) sign at 
the right hand side refers to a brane at the left (right) of 
the bulk interval. This equation does not provide any ad- 
ditional information with respect to the bulk equations. 
We can indeed rewrite the right hand side in terms of bulk 
quantities using eq. (|13p : the resulting equation is simply 
the constraint equation ((22|) at the brane location. The 
linearization of ^ is instead nontrivial: 



= ±^U,,S(b, , (26) 

branc Z 



□Q 



(n) 



(") 



(29) 



Eq. (|28l) is the decomposition in Kaluza-Klein modes; 
each mode is characterized by a wave function Wj-"^ in 
the bulk, as well as a 4D (quantum) field Q*^"^. Equa- 
tions (l24l) and (l27l). with the substitution v. 



Vi and 

□ — > m^j, provide the complete eigenvalue problem to 
determine the eigenmasses and the bulk profiles of the 
scalar modes of the system. It is clear from the form of 
the equations that this problem is well posed, and can be 
uniquely solved. We discuss this in detail in Section IVTl 
and provide an explicit example in Section IVlIl 



C. Boundary conditions for (infinitely) stiff brane 
potential 

One can obtain a simpler set of boundary conditions 
than (|27p in the limit of infinitely stiff brane potentials. 
Let us Taylor expand the brane potential for small fluc- 
tuations S4>i around the background values 4>i- 



where again the left hand side is evaluated in the bulk 
immediately next to the brane. 

The role of the boundary conditions is to complement 
the bulk equations ()24p and form an eigenvalue problem 
that can be immediately solved to obtain the physical 
scalar excitations of the system. By looking at eqs. 
we see that the most useful form in which these equations 
can be written is v[ = CijVj, where the coefficients Cij 
depend on background quantities. Considerable algebra 
is required to obtain this equation starting from 



U [if^ = 0,-1- 54>^) = U + U^^ 6(j)i + hj,,j54>^ 6(1) j + , 

(30) 

where the potential and its derivatives on the right hand 
side are evaluated at the background solution (j)i . As can 
be seen from and ([13]), only the expectation values 
of the brane potential U and its first derivatives U^i are 
relevant at the background level. The brane "masses" 
niij = U^ij only enter in the boundary conditions for the 
linear perturbations; higher-order terms in (I30p are in- 
stead relevant only beyond the linearized level, and can 
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be disregarded in our study. The stiff potential limit is 



the limit for which the "masses" 



are much greater 



than any other mass scale in the problem. In the original 
Goldberger-Wise [ll| stabilization mechanism, this is the 
limit of large A^^^; it is also explicitly noted there, that 
the equations considerably simplify in this limit. This 
limit can always be imposed by simply adding a quadratic 
potential term AJ7 = Ay {(pi — 4>i) {ipj — centered on 
the background solution, and then taking the limit of 
large Ay . Adding this term does not modify the back- 
ground solution. 

The simplification occurs because in this limit Scjii — > 
at the brane location; this can be seen from eqs. (1^51) . 
More precisely the fluctuations Scfii on the brane "ad- 
just themselves" to the O {U~^) level required to sat- 
isfy We actually do not need the explicit solutions 
for the scalar fluctuations. If we consider the equations 
in the original set of variables {5(j)i, $} we see that the 
boundary conditions (j26p are the only terms in which 



the second derivatives U^ij are present. These equations 
can be solved for sufficiently small brane, but then, 
once they are satisfied, the only role that these equations 
play in the stiff limit is to impose that <50i^branc at the 
boundary can be set to zero in all the other equations of 
the system. The situation is completely analogous to the 
scattering of light massive objects against an infinitely 
heavy object. An object of large mass A4 acquires an in- 
finitesimally small, 0{1/M), velocity in the scattering. 
For M — > oo we simply disregard the motion of the heavy 
object, and the value of M drops from the problem; the 
role of the heavy object is to ensure that the momentum 
conservation equation is satisfied, but then this momen- 
tum conservation equation plays no role in the dynamics 
of the light mass(es) participating in the scattering. Eqs. 
(PS)) are analogous to the momentum conservation equa- 
tions, U^ij is the analog of the mass M, and 8(j)i, 

ibranc are 

analogous to the velocity acquired by the heavy object. 
In this limit, eq. (j23p then imposes the condition 



Vi.hr. 



= N4>[ 



2, brane ' 



(31) 



where the proportionality constant A/" is the same for 
all modes. Eq. ([51]) provides a system of A — 1 inde- 
pendent boundary conditions in the eigenvalue problem. 
The reason why the number is A^ — 1 (rather than A^) is 
because the overall normalization of a mode - which is 
proportional to M - cannot be specified by the linearized 
problem we are solving (if one multiplies the solutions 
of a linear system by a common factor, one still has a 
solution). Fortunately, we do not need to know M if 
we are only interested in solving the linearized problem 
for the eigenmasses of the modes, and therefore we can 
simply set M to any convenient nonvanishing value (see 
Section Ivn. ^ 



The "missing" boundary condition, in addition to the 
A^ — 1 independent boundary conditions (1311) , is obtained 
by evaluating the constraint equation (|22|) at the brane 
location, for (5(/)i brane = : 



2A' 



= 



(32) 



We need to rewrite this equation in terms of v[ and M . 
Using (pij) and 50^, brane = in eqs. we obtain 



y/2A' 

■ ^5/2 



(33) 



To obtain the desired expression for $' Ibranc, we then dif- 
ferentiate with respect to z the second of the bulk equa- 
tions (| A-4p . ^ Using ([55]) in the resulting expression, we 
obtain a relation for $' Ibranc in terms of v[ brane ^"^"^ 
Using this relation, and eq. ((33)) . we can then cast the 
condition ([5^ into the desired form. In terms of the 
bulk wave functions, the condition reads: 



brane 



5 A' , 
'2 A' 



A' 
A 



A 



~^ 3M^A 



In the stiff brane limit the A^ equations ((3T|) and ([34] 
replace the conditions (^7)). 



D. Boundary conditions without a boundary brane 

Next we comment on the possibility that one or both 
boundary branes are absent. Assume that the bulk coor- 
dinate extends to infinity in that particular direction(s). 
In this case, the boundary conditions for the perturba- 
tions can be dictated by the specific problem under con- 
sideration. A typical requirement is the one of normal- 
izability. The set of differential equations ([24)1 has 2A^ 
solutions. In the example that we study in Section rVIIi 
one finds that half of these solutions exponentially grow 
at z — > oo, while the remaining half exponentially de- 
crease. Therefore, eliminating the exponentially growing 
solutions provides precisely A^ boundary conditions, as 
it was the case for the boundary conditions enforced by 
a boundary brane. Other models are characterized by a 
horizon at some bulk position, and one then typically re- 
quires that the solutions should be purely infalling modes 
at the horizon. Also this requirement corresponds to A 
boundary conditions. 



The normalization of the modes, or, equivalently, the value of 



Af in eqs. JSTJ, is needed if one instead wants to compute the 
couplings of the perturbations to each other, or to other fields 
(because the couplings are determined by the actual value of the 
wave function). The normalization is obtained by canonically 
normalizing the modes in the quadratic action of the perturba- 
tions, see Subsection IIII El 
^ Note that eq. l|33|l is a boundary condition, and so we cannot 
obtain an expression for "I"' Ibranc by simply differentiating it. 
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E. Quadratic scalar action, and normalization of 
the scalar modes 



To properly normalize the scalar modes, we compute 
the kinetic term of their quadratic action, obtained by 
expanding the starting action ^ at second order in these 
perturbations. We find 



(fx 



-ViUVi +dy \. 



AA' 



(35) 



If we replace ^' in the first line of this expression through 
the constraint equation (|2ip. ^t* = — ^, we can imme- 
diately see that the kinetic term is manifestly positive 
(recall that 77™ = —1), which ensures that the scalar 
system has no ghosts. The second line of (|35|) has been 
obtained following the steps outlined before eq. (19) of 
, where an analogous computation was performed for 
the case of a single scalar field. By decomposing Vi as 
in (UHl), and, analogously, $ = E„ (z)Q(") {x), we 
arrive at 



S, 



2, kin 



where 



Cmn =2 I dzv. 



(™) ,,(") 



3M3A4 
AA' 



([)(™) <J)(n) 



(36) 



(37) 



In evaluating the boundary term, one can make use of 
eq. (|A-7|) to express ^^"^ in terms of w-"'' and w-^"''. 

We note that the result ([57]) is the most immediate 
generalization to N fields of the expression (26) obtained 
in for the single field case (we note that the sign of 
the boundary term in the intermediate expression in eq. 
(26) of [T4| is incorrect). Hermiticity of S2 ensures that 
eigenmodes with different mass are orthogonal, so Cmn oc 
5mn- Imposing C„„ = i, we then recover a diagonal and 
canonically normalized kinetic term 



S- 



2,kii 



^ \ I d*a;Q(")nQ("). (38) 



If the linearized equations can be solved analytically, 
one can leave the normalization of the modes (specifi- 
cally, the quantity Af, in the case of a stiff boundary 
potential), and then fix it through C„„ — ^. li the equa- 
tions can be only integrated numerically, one needs to 
set a provisory value for the normalization by (arbitrar- 
ily) fixing the value of one of the wave functions at one 
boundary (for instance, in the example that we study in 
Section IVIll we set = 1 at the UV boundary); after 
performing the numerical integration, one can then in- 
sert the solutions in (1371) and obtain the provisory result 

= D^"'' provides 



The rescaling v. 



(») 



the correctly normalized modes. 



IV. VECTOR PERTURBATIONS 

For the vector modes, the 5/i and fii> components of 
the linearized Einstein equations ^ in the bulk give, 
respectively. 



□Bp = 



A' 







(39) 



while the 55 component trivially vanishes. Moreover, the 
linearization of equations (0]) has no contributions from 
the vector modes. The bulk equations therefore enforce 



A-'^ (z), with Dbf, = 0. This imme- 



diately implies that there are no Kaluza-Klein modes in 
the vector sector, except for a zero mode. 

If branes are present, it is also immediate to see that 
the linearization of the boundary conditions for the scalar 
fields, eqs. ([5]) have no contributions from the vector 
perturbations. The linearization of the Israel junction 
conditions, eqs. ©, is instead nontrivial. We see from 
([7]) and (|A-2p that the brane stress-energy tensor is 



B^— independent. Using (IA-3|) . we therefore have 



(40) 



Under the assumption of a Z2 symmetry, B^ needs to 
be odd across the brane,'' which then enforces = 0. 
The absence of vector modes agrees with the counting of 
the number of physical degrees of freedom given in the 
Introduction. 



V. TENSOR PERTURBATIONS 

The pLV components of the linearized Einstein equa- 
tions ([2]) in the bulk give 



h'' 



A' 

3^ K 

A ^' 







(41) 



while the other components trivially vanish. Moreover, 
the linearization of equations (|4]) has no contributions 
from the tensor modes. 

If branes are present, we again find that the lineariza- 
tion of the boundary conditions for the scalar fields, eqs. 
([S]) have no contributions from the tensor perturbations. 
The linearization of the Israel junction conditions, eqs. 
© instead gives 







(42) 



* The simplest way to see this is to consider the gauge invariant 
combination out of the two vector perturbations appearing in 
the metric II14I I. As shown in II17I I. — > _Ep — under the 
infinitesimal coordinate transformation IjlSj l. One can also show 
that, under this transformation, — >■ B^^ — Therefore, the 
gauge invariant combination is Bfj, — E'^. As E'^ is odd across 
the brane, this implies that is also odd. 
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We note that the tensor modes only "respond" to 
the background geometry, and not to the details of the 
sources. Therefore, these results are similar to those al- 
ready derived for the case of a single scalar in the bulk 
(see, for example, [23 )• F^^f. [201 the Kaluza-Klein 
eigenmasses were shown to be nonegative, so that there 
is no instability in the tensor sector. For completeness, 
we summarize this computation here: 

One starts by decomposing the 5D tensor field as 



where satisfies the Schroedinger-hke equation 



h'L 



3A" 3A'2 



One can further define the operators 



(43) 



(44) 



(45) 



so that the bulk and brane equations become, respec- 
tively, 



2?_ hr. 







branc 



, (46) 



where in the last step we have used the Z2 symmetry 
across the brane. 

If the first equation of p6|) is multiplied by /i* from the 
left and integrated over the bulk coordinate z, the left 
hand side can be then integrated by parts. The result- 
ing boundary term then vanishes because of the second 
equation of (|46)) and leads to the condition [20l| 



2 Jdz\V^ fc„P 
m„ ~ - — : T^-T — > 



J dz\h„ 



(47) 



Using Eqs. (|46p we can also immediately determine 
the existence of a tensor massless mode, characterized 
by the bulk profile 



(48) 



From (|43p and the metric decomposition (1141) . we re- 
cover the well-known fact that the massless tensor mode 
has an identical bulk profile as the background geometry. 
Incidentally there is no issue with using the decomposi- 
tion p4p for the massless tensor mode because the tensor 
equation of motion does not change JJj. 



^ Notice that [20l | defines X>± in the opposite way, since their f2 is 
l/A. 



VI. EIGENVALUE PROBLEM 

In the three previous Sections we have obtained the 
canonical modes in both the scalar and tensor sector, 
while we have shown that there are no physical vector 
modes. We have decomposed the canonical perturba- 
tions in a Kaluza-Klein sum, and obtained the equations 
satisfied by the wavefunction of the modes. In the two 
following Subsections we outline the eigenvalue problem 
that can be solved to obtain the properties of the physical 
modes. 



A. Scalar sector 

When iV bulk scalars are present, the scalar sector of 
the perturbations is characterized by the N physical 5D 
perturbations u^, defined in (j23p . The wave functions of 
the corresponding Kaluza-Klein modes satisfy the TV sec- 
ond order differential equations p4|. Each KK mode is 
therefore characterized by IN -I- 1 parameters (the mass 
m„, and the IN values required to specify the Cauchy 
problem), so that the bulk differential equations need to 
be supplemented by 2A^ 4- 1 conditions. Each bound- 
ary brane enforces iV conditions, given by eq. (|27p in 
the case of finite brane potentials, and by eqs. pip and 
([M]) in the case of infinitely stiff brane potentials. We 
discussed in Subsection IIII Dl the typical boundary con- 
ditions that can be imposed in the case that one or both 
branes are absent (in general, we expect N conditions 
per boundary). One additional condition is obtained by 
fixing the overall normalization of the modes. One can 
typically start by fixing a provisory (and, generally, in- 
correct) normalization; for instance, one can require that 
one of the wave functions evaluates to 1 at one bound- 
ary. This, together with the 2N conditions coming from 
the two boundaries, allows the system of linearized equa- 
tions to be completely solved. In this way, one obtains 
the masses of the physical modes, and the wavefunctions, 
up to an overall, yet to be specified, normalization. We 
stress that the overall normalization of the solutions is 
an irrelevant quantity in a linearized system of equations 
(if we rescale all the modes of a solution by a common 
factor, we still have a solution). Therefore, the system of 
equations can be solved, and gives the correct values of 
the eigenmasses, for any arbitrary overall normalization. 
Still, the overall normalization is crucial to determine the 
couplings of the perturbations, (since the couplings are 
determined by the values of the wavefunctions) ; Subsec- 
tion |IITE] explains how to rescale the solutions, so as to 
obtain the correct normalization, once the linearized sys- 
tem has been solved. 

In many cases, the bulk equations cannot be solved 
analytically; we expect this to be the norm in the scalar 
sector, where the perturbations are coupled in the equa- 
tions (this happens even if the bulk potential is a sum of 
separate terms, due to the mixing of the scalar field per- 
turbations with the metric perturbations). In this case. 
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the eigenvalue problem needs to be solved with a shoot- 
ing method; one fixes half of the 2N parameters that are 
necessary to determine the bulk evolution, by enforcing 
the N conditions at one of the two boundaries. Next 
one guesses the remaining N parameters, and solves the 
bulk equations which enable the wave functions to be 
evaluated on the other boundary. If the resulting so- 
lutions happen to satisfy the boundary conditions also 
there, then one has obtained a physical mode of the sys- 
tem. Typically, the initial guess is not correct, and one 
needs to employ some numerical scheme to obtain the so- 
lutions. For example, one can compute by how much the 
wavefunctions evaluated at the second brane differ from 
the expected boundary conditions, as a function of the 
initial guesses. An TV-dimensional Newton's method can 
then be implemented to find the zeros of this function. 
We perform this algorithm in the example studied in the 
next Section. 



B. Tensor sector 

The tensor sector is significantly simpler than the 
scalar sector. The wave functions defined in equa- 
tion ((43|) . satisfy the second order differential equation 
in the bulk. Each solution is in principle character- 
ized by three parameters: the mass to„ of the eigenmode, 
and two integration constants Cim C2n- However, as for 
the scalar case, only the ratio of these two constants, and 
not the overall normalization of the solution, can be de- 
termined from the linearized problem. Therefore to just 
obtain the cigenmasses and bulk profiles we can simply fix 
an arbitrary normalization by requiring that /i„ acquires 
a nonzero (but arbitrary) value at a given position (typ- 
ically, at one boundary brane). This gives one condition. 
The other two conditions are enforced by the boundary 
branes (each brane enforces one condition, given by the 
last expression in pS)) ). or if a brane is absent, by the re- 
quirement that the wavefunction is normalizable. These 
three conditions then allow to„,Ci„, and C2n to be de- 
termined with the bulk profiles known up to an overall 
normalization. The correct normalization can be then 
obtained from the quadratic action of the tensor modes, 
analogously to what we dicussed for the scalar sector. 



A. The 5D model 

We review the soft wall background solution ^ISj (with 
the only difference that we use our convention for M'^, 
that corresponds to in [l^). The model is charac- 

terized by the two fields (j) and T, with the bulk potential 



1 



1 



V = -zy^rT^e6^^"(i+-) + -v^k 



(1 + ^) 



V2 



. g 12MJ(1 + i^) 



M3/2 



V V2 \/6 
One obtains the background solution 



. V3 Af 3/2 



(49) 



e 3 



k z 



A{z) 

T{z) = -2^2^1 + l/vM'-'''^ {jiz) 



(50) 



The parameter is a dimensionless constant, while fj, is 
the soft-wall mass scale. As the bulk volume diverges 
at 2; — 0, a UV brane is placed at zq = 1/k. The po- 
tential U on this brane is chosen so that the solutions 
([50| satisfy the boundary conditions there: specifically, 
the background values of U and U^i are determined from 
and p^ . We then assume that the brane potential 
contains large quadratic terms (j30p . so that the bound- 
ary conditions for the perturbations can be given in the 
infinitely stiff limit of Subsection IIII CI 

The scalar field configuration provides a finite bulk ge- 
ometry in the limit z — 00, so that there is no need to 
include a boundary IR brane at large z (in fact this is 
the reason for why it is a "soft wall", as opposed to a 
sharp brane or hard- wall cut-off). As we will see in the 
next two Subsections, the requirement that the perturba- 
tions are normalizable there provides sufficient boundary 
conditions to fully determine them. We focus our study 
of the perturbations of this model to the choice ly = 2, 
as this gives rise to a linear Regge-like mass spectrum, 
TOjj oc n, for the KK modes, which is similar to that 
encountered in the hadron spectrum of QCD. 



VII. AN EXPLICIT EXAMPLE: THE 
DYNAMICAL SOFT WALL 

As an explicit example of the general method that we 
have outlined above, we will consider the dynamical soft- 
wall solution found in Ref. [l^. In addition to the met- 
ric, this solution involves two bulk scalar fields. It pro- 
vides a dynamical realization of the holographic soft-wall 
model for QCD [2l|, [g^], as well as applications to the 
electroweak sector of the Standard Model [23[ . 



B. Tensor modes 

The tensor modes for the model were already studied in 
[2^ . We summarize these results here for completeness, 
and to provide an example of an eigenvalue problem that 
can be solved analytically (and that is technically simpler 
than the problem for the scalar modes studied in the next 
Subsection) . 

It is convenient to use the variable /i^^ = JgJJ/A^ {z) 
(where TT denotes the transverse-traceless component) 
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for the eigenvalue problem. The Kaluza-Klein decompo- 
sition is analogous to , and we denote by /i„ the wave 
function of the n— th mode (/i„ = A^"^^^ hn, where ft,„ is 
introduced in (|43|)). In the background (|50|) . the bulk 
equation eq. (HIT) becomes 



(Phn 

dP 



- 4z 



'i]^+fhlh^^O, (51) 
z J dz 



where we have introduced the dimensionless quantities 
z = pLZ and m = m//i. The boundary condition at the 
UV brane then has the form 



dhn 
dz 







Z—fl 



(52) 



(where jl = fx/k corresponds to the location of the UV 
brane) while the normalizability requirement at z — > cx) 
translates into the requirement that the solution de- 
creases sufficiently fast at large z. More precisely, since 
the wave function of the canonically normalized mode is 
hn = A^/"^ hn, we require that J dzA^ h\ < oo. 

The normalizable solution of Eq. (|5T|) is (up to a nor- 
malization constant) the Kummer's confluent hypergeo- 
metric function 



U 



-1, 2z' 



The eigenmasses are obtained by imposing (I52p . In the 
limit /i = /(i/fc<Cl, we can expand ((55)) at small z, and 
obtain 



const. 



16z 



;2 



(-8 + m2)r(-5a) 



+ 0{i^) . (54) 



Eq. (15^ is approximately satisfied at the poles of the 
gamma function, namely for 



2V2n 



0, 2, 3, 



(55) 



(we verified numerically that there is indeed no physical 
mode corresponding to n = 1). Wc note that m^^ > in 
agreement with the general result (|T7|) . 



Scalar modes 



For the background solution ((50| . and in terms of the 
dimensionless quantities z, m„ defined in the previous 
Subsection, the bulk equations (1^^ become 



-9 + [-9 + 16z2 (-3 + ^2 + 2^4)] 



452 (3-f 4i2)2 



16z(-9-f 24z2-M6z4'> 



-V2 



v''=l -m2 



3(3 + 4z2)" 
243 + 8z2 [81 + 8z2 (l71 -I- 78f2 + 4^4)] 



36z2 (3 + 452)" 
16z (-9 -I- 2452 + 16^4) 



U2 



3 (3 -1- 4z2 



-Vi 



(56) 



where prime (') denotes differentiation with respect to z, 
and where the subscript 1 (2) corresponds to the field (j) 
(T). As we discussed in Section ITTll to fully determine a 
mode we must specify the five quantities (2A^ -I- 1, with 
N = 2) 



Vl\ 



V2\ 



(57) 



where we recall that jl is the position of the UV brane in 
the rescaled variable z. 

We assume a stiff brane potential on the UV brane, 
which enforces the Dirichlet boundary conditions Scf) = 
6T — (see Subsection IIII C|) . This in turn results in the 
two boundary conditions ((3T|) and (|34l) . Moreover, we fix 
the arbitrariness of the overall normalization by fixing 
the value of W2 = 1 at the UV brane. This, together with 
Eq. (ini) gives 



^^1 



■3^' 



V2\ 



(58) 



We are left with the three parameters w^l/i. 



(53) subject to the constraint 



:^-2Ia = "0(3 + 4m^ 



2/i ""'^ 
1 



up. 



:^ (9 + 16 /i*) + 4 



3 + 4/i- 



r/2 



(59) 



We used this constraint to determine v'i\f^ as a function 
of V2\fi. and ni^. 

The bulk coordinate z extends to -f 00. To understand 
the role of the associated boundary conditions, we stud- 
ied the bulk equations in the limit of large z: 



■j'l = [4/z2 -f O (z")] t-i 
16^3 



16/i^ 



z + 0(z 



-3^ 



V2 , 



-z + 0(z-3) 



Vl 



4^4 



z2 + 0(z") 



V2 



(60) 



We obtained the approximate solutions to these equa- 
tions under the assumption that one mode is significantly 
larger than the other one in this limit; the subdominant 
mode can then be disregarded in the equation of the dom- 
inant mode; this equation can be solved analytically, and 
we can then insert this solution in the remaining equa- 
tion to obtain the subdominant mode. We then studied 
the large z limit of the solutions, and verified that the 
starting assumption (the subdominant mode can be ne- 
glected in the equation of motion of the dominant one) 
holds. In this way, we obtained four solutions, that form 
a complete basis for the solutions of the bulk equations. 
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At z ^ 1, the four solutions read 



V2 



V2 



Vi 



Vi 



3Ci 
3Di 



e < ui ~ ^ e 
Vz 

e"^ < wi ~ ^ e^ , 



C2 

3-D2 iy2 i;2 

^^3 «-2-^e3 . 



(61) 



We note that two of these solutions are exponentially de- 
creasing at large z, while the other two are exponentially 
increasing. 

As we discussed in the previous Section, at this stage 
of the computation we can only guess some values for 



the initial parameters W2|/i, 



This guarantees that 



the boundary conditions at the UV brane are satisfied. 
If we start from any point in this two dimensional space 
of initially guessed parameters, and solve the bulk equa- 
tions ([55)) from the UV brane to the asymptotic z ^ 1 
region, we obtain solutions whose large z asymptotics is 
a linear combination of the four modes (I6ip . Only lin- 
ear combinations that have Di = D2 = correspond to 
a normalizable, and hence physical, Kaluza-Klein mode. 
Given that we need to satisfy two conditions, and the 
space of initial parameters has dimension two, we expect 
that only a discrete set of points in this space corresponds 
to a physical mode. Our goal is to identify these points. 

As can be expected, none of the initial guesses cor- 
responds to a physical mode, and therefore we resorted 
to the following numerical algorithm. For the present 
discussion, let us denote by a = and /3 = V2\ji, the 
two coordinates in the space of initial parameters, and by 
/ (a, P) and g (a, /3) the values of the two wavefunctions 
Vi and V2, respectively, at some large value Zend- Physi- 
cal modes correspond to zeros of these two functions. We 
numerically searched for zeros using Newton's method: 
starting from an initial guess {olq, /3o} we initiate a suc- 
cession {a„, /3m} that, if it is convergent, converges to a 
zero of / and g. The iteration step of the succession is 



ttn+l 



/,/3 .9 - / 

f,a g,i3 - f,p g,t 
g,a f - f,ag 
f,a g,i3 - 1,0 g,c 



(62) 



where the subscript comma denotes differentiation us- 
ing finite differences, namely = ■^("+'^'^)~-^("~'^'^) ^ 
and similarly for the other derivatives. The iteration 
step is obtained by Taylor expanding / (ctn+i, Pn+i) and 
g(oin+i^ Pn+i) in terms of /, g and their first deriva- 
tives, evaluated at {a„, /?„}. The expressions (|62l) are 
the algebraic solutions of these equations, after setting 

/ (a„+l, Pn+l) = g (an+l, Pn+l) = 0. 

Ideally, the algorithm converges to real solutions only 
for Zend = oo. In practice however, the numerical prob- 
lem we are solving is rather challenging. Indeed, we need 



to identify exponentially decreasing solutions among ex- 
ponentially growing ones. No numerical solution starts 
from values of {a, /3} that exactly correspond to a solu- 
tion. The discrepancy, however small, necessarily results 
in an exponentially growing solution at z 3> 1. The larger 
the value of Zond , the smaller the initial discrepancy needs 
to be, if one hopes that the growing mode is still subdom- 
inant at Zend- After some trials, we found that values of 
Zend — 5 result in successions that converge after ~ 10^ 
or 10^ iteration steps of ([62]) . For larger values of Zcnd, 
the successions typically keep spanning large areas of the 
{a, j3} plane without showing any sign of convergence. 
We note that Zend — 5 are the values for which the ex- 
ponential suppression/enhancement in (jFl]) starts to be 
significant. 

We generally expect that the effect of having a finite 
value for Zend, rather than Zond = oo, is the following: any 
zero of / and g that we find arises because the exponen- 
tially decreasing and the exponentially growing modes of 
(j6ip provide an equal contribution to the wave function 
at z ~ Zend (if tbis was not the case, the wave func- 
tion would exponentially decrease, and not vanish, at 
^cnd)- For sufficiently large Zend, this however implies 
that the growing mode is extremely subdominant at the 
UV brane, and that the initial choice of W2|/i, is very 
close to the physical one. For instance, in obtaining the 
eigenmasses that we show in Figure [5l we tried to push 
Zend to the highest possible values that allowed for a con- 
vergence of ([62]) , and we verified that a small decrease of 
Zend from this highest possible value did not significantly 
change the values of ^21^, to which ([62]) converged 
to. 

Our main goal in the current example is to study 
whether the soft wall model |15| is stable. Clearly, our 
numerical search will never produce a definite proof of 
this, since one may in principle think that modes with 
negative m? exist, but our numerical scheme could not 
find them. However, as we now discuss, our search al- 
gorithm converged to more than 100 eigenmodes with 
TO^ > 0, and there were no modes with < 0. We 
believe that this provides substantial evidence that the 
soft wall model is stable. 

We performed two extensive searches, with different 
values of //^. In the first search, we fixed ji — 0.01, and 
^cnd — 5, and we took a grid of 250 x 500 values in 
the plane of initial values; specifically, the values of to^ 
ranged from —49.9 to —0.1 in steps of 0.2, and W2|/i from 
—49.9 to 49.9 in steps of 0.2 (starting from negative val- 
ues of so as to maximize the chance of finding some 
tachyonic solution, if it existed). Starting from each of 
the 125, 000 points in the grid, we performed 400 itera- 
tions of ([5^ . The iterations converged to 118 final points 
with TO^ > 0, and with \/ P -\- g"^ <C 1 (therefore, they 
correspond to stable solutions of the problem). In the 
second search, we fixed fi = 0.001, and Zend = 5- We 
took a grid of 100 x 200 values in the plane of initial 
values; specifically, the values of ranged from —59.7 
to —0.3 in steps of 0.6, and Ual/i from —597 to 11343 in 
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steps of 60. We performed 400 iterations of ((62)) . which 
converged to 67 final points, with > 0, and with 

^« 1. 
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FIG. 1: Values of = m^//i^ and U2,o = ^'21/5 to which the 
algorithm described in the main text converged (with fi = 
0.01). Note that all the solutions have positive mass-squared, 
m^. Also shown is the linear fit along which the parameters 
of the solutions are approximately distributed. 

The values of and Wjl/i of the convergence points 
are shown in Figure [U for the search with /2 = 0.01. We 
see from the Figure that many of the solutions have a 
final value of these parameters outside the range of the 
starting grid. This shows that the algorithm is able to 
find solutions in a broad region of the parameters. We 
see that the values shown in the figure are roughly dis- 
tributed along the line, iijIa — 25m^ — 42, although some 
scatter is present (for the search with jl = 0.001, we 
found that the parameters of the solution are roughly 
distributed along the line — 250m^ — 520). We also 
notice that the solutions we have found are not equally 
spaced along this line: this is a clear sign that the search 
has not found all the modes with the parameters in the 
range shown in the Figure. Particularly, it is natural to 
expect that the method failed to obtain modes with val- 
ues of and V2\ji very different from those in the initial 
grid. While the exact spectrum of the soft wall model [lH| 
is by itself a very interesting subject worth studying, this 
is not the goal of the present work. The main purpose 
of Figure [1] is to show that the algorithm that we have 
implemented did not find any unstable (m^ < 0) mode. 

In Figures [5] and [3] we show the bulk profile of vi 
and V2 , respectively, of the three lightest modes obtained 
from this search. These wavefunctions are shown with 
the "provisory" normalization set by eq. ([55]) . In Fig- 
ure m we show instead the bulk profile \/ v\ + V2 of the 
three modes after they have been properly normalized 
(see Subsection nil E[) . As typical, the profile of modes of 
increasing mass presents an increasing number of maxima 
and minima, and extends further in the bulk. 

The hierarchy in the model (namely, the ratio between 




FIG. 2: Bulk wavefunctions of v\ for the three lightest scalar 
Kaluza-Klein modes for fi = 0.01. The corresponding eigen- 
masses are ~ 0.41, 3.47, 6.25. 
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FIG. 3: Bulk wavefunctions of V2 for the same modes appear- 
ing in Figure [21 



the IR and the UV mass scale) is given by Jl. While a 
solution to the hierarchy problem requires jl ~ 10~^^, the 
two searches we have performed are limited to very mod- 
erate hierarchies, namely /2 — 0.01 and 0.001. However, 
we also studied how the lightest eigenmasses scale with 
/2, and we believe that our results can be extrapolated to 
small hierarchies. In Figure [5] we present the value of four 
light eigenmasses for jl as small as 10~^. We find that 
the lightest eigenmode shown behaves differently than 
the heavier ones. For all these other modes, the ratio 
TO^ ~ m^l approaches a constant value at small . 
Instead, we find that the numerical values obtained for 
the lightest eigenmass are extremely well approximated 

by 



^radi( 



1.9 



(63) 
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FIG. 4; Properly normalized bulk profile of the same modes 
appearing in Figures [5] and [31 
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FIG. 5: Values of the lightest eigenmasses as a function of 
fi = fi/k. 



those obtained in the Goldberger-Wise llll s tabilization 
mechanism, as shown in the analysis of [12|. [2^ (in the 
computation of 112], the logarithmic suppression is en- 
coded in the factor u/k^l/ (kro) ~ 1/log (UV/IR)). 

The nature of these modes as IR modes is also con- 
firmed by the computation of their normalization, that 
we performed for the specific case of /i = 0.01. Specifi- 
cally, we evaluated C„„ in eq. (I57|) for the four modes of 
Figure [5] We found that the contribution to C„„ from 
the bulk integral is about 5 orders of magnitude greater 
than the contribution from the UV boundary at z — fl. 
Since the wavefunctions Vi (and, thus, their bulk inte- 
gral) are peaked at values of z of the order of the inverse 
IR scale, this confirms the IR nature of the modes. In an 
analytic computation, one would expect no contribution 
to C„n from the boundary at z = cx). In our numerical 
computation, the contribution from the IR "fictitious" 
boundary at Zcnd = 5 ranges from being 5 to 7 orders 
of magnitude smaller than the bulk integral. This is a 
measure of the goodness of our numerical results. 

A second check on our numerics is provided by the 
orthogonality of the modes. Ideally, the modes should be 
perfectly orthogonal to each other: Cmn = for m ^ n. 
We find that s^n = Cmnl \/Cram Cnn ioT m ^ Ti is of 
O (lO^^) for the four modes of Figure[S] As a final check 
of the numerics, we verified that the wave functions Vi of 
the four modes of Figure [5] is in good agreement with the 
analytic large z solution that we wrote in the third line of 
eq. (|5T|) . We recall that normalizable modes must behave 
as linear combinations of the first and third line of (j6ip 
at large z. Since the mode in the third line (cx e~^ 
decreases less than the mode in the first line (oc e~^ ), 
only the e~^ component will be visible if the mode is a 
generic linear combination of these two solutions at large 
z. 



VIII. CONCLUSION 



We were actually able to obtain this eigenmass for values 
of jj, as small as 10"^" (for smaller values, the algorithm 
was no longer able to converge to the solution), and we 
verified that the fit (|63| continues to be valid also at such 
low jl. 

In eq. (1631) we referred to this mode as a radion, be- 
cause the behavior of its mass is consistent with that 
of the radion in the Goldberger-Wise [llj stabilization 
mechanism, while the behavior of the heavier masses is 
consistent with the behavior of the higher KK modes in 
[T]| . We recall that /i is the IR scale of the model, while 
k is associated with the UV scale. We therefore find 
that the wave function of the radion and of the higher 
KK modes are peaked at values of the bulk coordinate z 
of the order of the inverse IR scale. Moreover, we find 
that the higher KK modes obtain a mass of order the 
IR scale, while the radion mass is suppressed with re- 
spect to this value by the logarithm of the ratio between 
the UV and the IR scale. These properties coincide with 



The main goal, and result of this work is to provide a 
formalism to study the perturbations, and the stability 
of codimension one brane compactifications with N bulk 
scalar fields. We obtained the closed set of equations 
for such perturbations, and arranged them in an explicit 
eigenvalue problem. These equations are valid for arbi- 
trary bulk and brane potentials for the scalar fields, and 
therefore can be readily employed to study any such con- 
figuration. Specifically, we identified the canonical per- 
turbations of the system, eqs. (1^51) . They satisfy the sec- 
ond order bulk equations (j24p . Such equations must be 
supplemented by boundary conditions. In Section IIIIBl 
and in the following two Subsections, we provided such 
conditions for the case of generic brane potentials, in- 
finitely stiff brane potentials, and no brane, respectively. 
As an example, we studied the system of perturbations in 
the soft- wall model of [l5| , which is characterized by two 
bulk scalar fields. We found the scalar perturbations in 
this model behave identically to those of the Goldberger- 
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Wise stabilization mechanism [ll| ; namely, there is tower 
of KK modes with mass at the IR scale, and a lighter 
radion mode, whose mass is suppressed with respect to 
the KK tower by a large logarithm; all these modes are 
peaked at IR values of the bulk coordinate. Note that 
our formalism assumes that there are no massless scalar 
modes, but can be straightforwardly generalized to study 
examples in this class. 
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APPENDIX A: Brane boundary conditions 

We give the definitions and intermediate steps relevant 
for the computation of the brane junction conditions, 
mostly following the notation of [23]. The brane posi- 
tion in the bulk is characterized by the vector {y^)i 
denoting the bulk coordinate of a point identified by 
the position on the brane. This defines the basis vec- 
tor e^j = ■ The normal of the brane is then defined 
by the requirement that it is orthogonal to the basis vec- 
tors, and normalized to one: ej^-j ua = and uati^ — 1- 
The induced metric and the extrinsic curvature if^i, 
are then 

It^i' = ^t^) ef^) [9AB - riA ub] , K^,,, = ef^) ef^j V^ns • 

(A-1) 

We are interested in the results up to first order in 
the perturbations. We use the metric (fT4|) . in the E = 
B = ~ gauge. We denote the brane position in the 
extra space as Zbrano = ^background position + C{x'^), whcrc 

is a perturbation. In vector notation, using the first 
component for the ordinary coordinates {A — 0, ... ,3), 
and the final component for the extra coordinate (A — 5), 



we have 
This gives 



{S^, C,^) and nA = A S'^, 1 + 



= [(1 + 2^-) 7]^^ + /ip 



(A-2) 



and 



A 



A'h., 



(A-3) 



APPENDIX B: Scalar equations in the bulk 

We linearize the bulk equations ^ and (HI in the scalar 
perturbations (fT4l) . in the E = B = gauge. As men- 
tioned in the main text, the off diagonal fj,iy Einstein 
equations enforce = — $/2, which we use to elimi- 
nate 5*. Among the remaining equations, the /j,5 Einstein 



equation is already given in eq. (|22l) of the main text. One 
can then show that the part of the fiv equation propor- 
tional to r/py (the one with the ellipsis in eq. ([201) ') can 
be written as a combination of (j22p and its z derivative. 
We can therefore disregard it. 

We are left with eq. (U) and the 55 component of 
which are, respectively. 



3A' 



-f30^$' + 2^2 + A%,j6(l)j = 



3M3 
2 

~3M3 



(A-4) 



where we have used ((22)) to simpliiy the second of these 
equations. The left hand side of these equations (jA-4|) 
will be denoted as Eq, and Eqgg, respectively and also 
the left hand side of eqs. ([^ are denoted as Eq,;. Using 
(P^ repeatedly, one can show that 



1 d A! 
2d;^^^^ + T^^^^^ 

-V2A3/2fEq 



1 



2A' 



' 3M3 
Eq55 



'fcEqfe 
Eq,; 







. (A-5) 



The first relation indicates that one equation in 
the {Eq, = 0, Eqgg = 0} system is redundant. The 
first and second relations together show that the 
{Eqj = 0, Eqjg = 0} system of equations for the + 1 
variables is equivalent to the system of equa- 

tions |Eqj = o| for the N variables {wi}. This is pos- 
sible because the variables {8(i>i^ $} are subject to the 
constraint (1^^ . 



APPENDIX C: Scalar boundary conditions 

In this Appendix we summarize the algebra needed to 
rewrite the boundary condition (|26p into the expression 
(P7)) . We recall that our goal is to write the boundary 
conditions as expressions of the form v\ = CijVj. 

We start by taking the z derivative of the definition 
of Vi, eq. ([23l) . We then eliminate S(j)' and from the 
resulting expression using ([^ and (f^ . Finally, we elim- 
inate S(j)i using (j23l) . This gives rise to the expressions. 



2A' {iM^A' 2 



± 



V2A' 2V2A 



(A-6) 



which are equivalent to the boundary conditions 
once we assume that the modes satisfy the bulk equations 
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in the bulk (which clearly must be the case). We need to 
eliminate $ in favor of a linear combination of w^. From 
the second equation of (|A-4|) , together with (f22|) and (f23l) , 
we find 



a linear combination of Vi only. This can be obtained if 



there is a second relation between 



is readily achieved if we multiply (|A-6[) by 
over i: 



and Vi. This 
and sum 



V2 



5 A' 



A 



'kVk 



A' 3M3 



(A-7) 



In principle, we can solve and obtain $ from this relation, 
and insert it in (jA-6[) . The resulting equations would 
be functions of Vi and v[ only, but would not be of the 
required form (the linear combination cx v'f. on the 
right hand side of (jA-7[) , means that a nontrivial inversion 
is needed to solve for v[). 

Instead to achieve our goal, $ should be expressed as 



iA' 
2T 



MPA' 

A'^'^'kVk 



''k<l>'k 



± 



y/2A' 2y/2A 



(A-8) 



Substituting this expression into (jA-71) gives rise to an 
equation containing only the modes $ and Vi. This equa- 
tion is then used to express $ as a linear combination 
of Vi. Inserting this expression into (|A-6p provides the 
boundary conditions (j27p in the required form. 
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